Evaluation of Sylvester type determinants using 
block-triangularization 

Olga Holtz * 

Abstract 

It is shown that the values of Sylvester type determinants for various orthogonal polynomials consid- 
ered by Askey in [J can be ascertained inductively using simple block-triangularization schemes. 

1 Introduction 

Richard Askey in [I] shows two ways, one matrix-theoretic and another based on orthogonal polynomials, 
to evaluate determinants 
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N x 2 () 
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D N+1 (x) := 







2 x N 
1 x 



which were first considered by Sylvester [5]. In addition, he obtains several generalizations of Sylvester's 
determinants and explores their connection to orthogonal polynomials. 

The purpose of this note is to show how the determinants from yQ can be evaluated in yet another way, 
based on partial information about left or right eigenvectors of the corresponding matrices coupled with 
a simple similarity trick. In all cases except one, only one (the most obvious) eigenvector is used. The 
exceptional case is the Sylvester determinant itself, where two eigenvectors are readily available and hence 
used to derive the result. 



2 Sylvester's determinant and two close variants 

Let us start with the Sylvester determinant. We want to prove that 



N 



Dn+i(x) = Y[(x + N -2j) 

3=0 



(1) 



(which is formulas (2.3), (2.4) from PP). Since the values of D\ and D2 agree with JTJ, it is enough to show 
that 

D N+1 (x) = (x-N)(x + N)D N _ 1 (x) (2) 
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for N > I. The Sylvester determinant -Dat+i is the characteristic polynomial of the matrix —Dn+i, i.e., 
D N+ i(x) = U\e*(D N+1 )( x + A )> where 
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Note that (1,1,1,....!) is a left eigenvector of Pjv+i corresponding to eigenvalue N and (1, — !,!,...,(— 1) ) 



is its left eigenvector corresponding to eigenvalue —N. The similarity transformation T>n^ 
where 



Tt 



N+l ■= 
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therefore reduces 2?jv+i to a block lower-triangular form 

" diag(iV,-7V) 

* Mn-i 



where 



M N - X := 



3 + (-N + l) -N + l 

N- 2 4 
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Now it remains to show that A^at_i is similar to Dn-x- But this is indeed so, since <Sjv-i-^ n-i<Sn-i = T^n-i 
where 
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This proves @ and therefore (Q. 



A related determinant 
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can be evaluated analogously and even more simply. We use the notation 
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for the matrix that satisfies Bn + \(x) = Y[\£a(B N+1 )( x + ■V)' 
To prove that 

jV 

B N +i(x) =]J[x+(N- 2j)a -N + j] 

(formula (2.8) from Q]), we will show that 

B N+1 {x) = (x + Na- N)B N (x - a). 



(3) 



(4) 



The vector (1, 1, 1, . . . , 1) is a left eigenvector of Bn+i corresponding to eigenvalue Na — N. The similarity 
transformation Bn+i > Tn+iBn+iT^+i where 



jv+i •= 



1111 
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10 
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reduces Bn+i to a block lower triangular form 

Na - N 





M N 



where 



-l-JV(a-l) 2a-N(a-l) 



-N(a-1) 

(N-l)(a-l) -2 3a 

(JV-2)(a-l) -3 

(JV-3)(a-l) 
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Finally, M. n is similar to £>at — a/: 



where 



A4at5at = S N (B N - al), 



'AT 











3 



This proves ©. 



As is shown in pQ, the formula J3J for B^ + i implies the formula JIJ for the Sylvester determinant D^+i, 
since 

n / x r B N+1 (ax) 

a — >oo d 

Also, (0) gives a formula for another related determinant 
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via the relation 



A N+1 (x) = 2 N+1 B N+1 (| 



3 Determinants for Krawtchouk and dual Hahn polynomials 

The determinant for the Krawtchouk polynomial Kn+i{x;p, N) is given by 



K N+1 (x;p,N) := 

-x+pN pN 

-z+piV + (l-2p) 

2(1 -p) 





p(iV-l) 
-a;+p7V + 2(l - 2p) p(N - 2) 



iV(l-p) -a;+p7V + iV(l-2p) 

With the usual notation (a)k for the shifted factorial (see ^1 formula (3.10)]), we need to prove (3.25)]: 



K N+1 (x;p, N) = (-x) N+1 . 

For that, it is enough to show that 

K N+1 {x- P ,N) = (-x)K N (x;p,N-l). 
As usual, we need to find an 'obvious' eigenvector of the matrix 
pN pN 

(1-p) pN+(l-2p) p(N-l) 
K. N+1 := 2(1 -p) pA^ + 2(l-2p) p(N - 2) ••• 



(5) 



(6) 



iV(l-p) pN + N{l-2p) 
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corresponding to eigenvalue 0. In this case, it happens to be the right eigenvector (1, — 1, 1, — 1, . . . , (— 1) N ) T . 
We then use the transformation 



N+l 





1 
1 
1 



* 
M N 



where 



pN + PN+(l-2p) p(JV-l) 

-pN + 2(1- p) P N + 2(1 -2p) p(N-2) 

pN 3(1 -p) 

-pN 



(-l) N pN 
Now, A4n is similar to the matrix ICn + I, viz. 

S N M N = (K N H 

via the transformation 

"10 
110 
110 
S N-= 11 



N 



(-1) N ••• 1 
to transform K-n+i into a block upper-triangular form 



(7) 







This proves © and hence (JSJ. 



N(l-p) pN + N(l-2p) 



(8) 



The determinant for dual Hahn polynomials, 



J Rjv + i(A(a;);7,5,7V):= 
X(x) + AT (7 + 1) 

(N + S) 
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where 

X(x) := -x(x + j + 5 + 1), 

appears in pQ as the first example beyond Krawtchouk polynomials in terms of simplicity of recurrence 
coefficients. This determinant seems much harder to evaluate, but this in fact requires just an additional 
shift of parameters. The formula for Rn+i is QJ (4.5)]: 

R N +i(K^);i,S,N) = (-x) N (x + j + S + l) N+1 . (9) 

It can be proved from the relation 

R N+1 (\(x)]<y, 5, N) = X(x)R N (X(x) + (7 + S + 2); 7 + 1, 5 + 1, N - 1). (10) 

The latter can be checked exactly as above, i.e., by applying the transformation 7/v+i of the form J7J) to the 
matrix 



1Zn+i(7,5) := 
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obtaining a block upper-triangular matrix 

^jv+i^iv+i(7! 5)T N+1 

where 
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and then transforming .Mjv(7, 6) m t° T^-n{i + 1, 5 + 1) + (7 + 5 + 2)1 by 

M N ('y,S)»S N M I f('y,S)Sji 1 , 
where Sat has the form (JSJ. This proves and 10. 
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4 Determinants for Hahn, Racah, and g-Racah polynomials 



The last examples appearing in yQ are three types of orthogonal polynomials, whose recurrence coefficients 
are fractional rather than polynomial. Recall that the recurrence coefficients appear as entries on the main 
diagonal and the first sub- and super- diagonals of the tridiagonal determinant that gives the value of the 
corresponding polynomial. The alternating sums of these coefficients are still zero, just like in the two 
examples we just considered in section 3. In other words, the determinants we now consider all have the 
form 
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: det(A(x) + Gn+i 



(11) 



where 

The coefficients a n and c„ for Hahn polynomials are 

(n + a + (3+l)(n + a+l)(N-n) 
a " ~ (2n + a + (3 + l)(2n + a + (3 + 2) 

n(n + a + (3 + N + l)(n + (3) 
Cn ~ (2n + a + f3)(2n + a + [3 + 1)' 

with X(x) = —x, while the coefficients for Racah polynomials are 

{n + a + l)(n + a + (3+ l)(n + 7 + 1)(N - n) 
Un ~ (2n + a + (3 + l)(2n + a + [3 + 2) 

—n(n + a + {3 + N + \){n + a + (3- j)(n + f3) 
Cn ~ {2n + a + P)(2n + a + P+l) 

with (3 + 7 + 1 = — N and X(x) — —x(x +7 + S), and the coefficients for q- Racah polynomials are 

(1 - afeg" +1 )(l - g" +1 )(l - q n - N ){l - cq n+1 ) 



(12) 



eg 



(1 -a6g 2 »+ 1 )(l -abq 2n + 2 ) 



(l-a6g 2 ")(l-a6g 2 " +1 ) 



(13) 
(14) 



with bdq = q~ N and X(x) = -(-q- x ){l - q x+1 cd) (see P Sec. 4]). 

Recall that Hahn and Racah polynomials are just limiting cases of g-Racah polynomials ^ Sec. 4], Pre- 
cisely, if we denote Hahn polynomials by _ffjv + i(A(x); a, /3, N), Racah polynomials by RApj + i(X(x); a, /3, 7, N), 
and g-Racah polynomials by QRA N+1 (X(x); g, a, 6, c, N), then we see that 



RA N+1 (X(x);a,/3,i,N) 

Hn+i{X{x); a, (3, N) = lim 



,. , 1 ^+1 QRA N+1 (X(x);q, q a ,qP, q~< , N) 
W" 1 ) (1 - qfN + 2 • 

RA N+1 {X{x); a, 0,j,N) 



7 — >oo 7 1 

So, it is enough to evaluate g-Racah polynomials. The formula is given in pQ in the form 

QRA N+1 (X(x);q,a,b,c,N) = {-l) N + l {q-*;q) N+1 {q x + l cd;q) N+1 , 
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where ( ; )k is the g-analogue of the shifted factorial ^ formula (4.13)]: 

fc-i 

3=0 

Its equivalent form that is more suitable for an inductive proof is 

N 

QRA N+1 (X(x);q, a, b, c, N) = J] (A(z) - A(n)). 

n=0 

This formula will be proved once we show that 

QRA N+1 (X(x); q, a, b, c, N) = X(x)q~ N QRA N (X(x -l);q, aq, b, cq, N-l), 

where A(a;) corresponds to the parameters (q, aq, b, cq, N — 1) so that 

cq 



\(x-l) = -(l-q- x+ ')(l 



bqN -f ) = q I A(X) + 1 + bqN 



1 cq 
q^bq^ 



(15) 



(16) 



(17) 



To prove let us start with an observation about our ansatz matrices Qn+i satisfying (jl 1|> and i|12|) . 
Suppose that such a Qn+i is transformed using the matrix Tjy+i given in J7J). Then, as we already saw, 
7/v+i^jv+i7/v+i is block upper triangular: 

" * 
M N 

Next, the transformation Sn given by (JBJ) reduces A4n to the tridiagonal form 



ao + ci a\ 

ci ai + c 2 a 2 
c 2 a 2 + c 3 a 3 



cat_ 2 ajM-2 + cjv-i 

Cat_i 







ajsr-i 

ClN-l + cn 



(18) 



Proving therefore amounts to showing that the matrix (|18H . with the entries a n and c„ coming 
from the determinant QRA N+1 (X(x); q, a, b, c, N), is similar to the matrix ^(Gn + A(x — 1)In) — ^(x)In, 
where Qn is an iVx TV-matrix of the form (| 1 1 1) — (| 1 . with the entries a n and c„ coming from the determinant 
QRA N (X(x — 1); q, aq, b, cq, N—l). (Note that, due to l(T7|l . the difference X(x) — -X(x — 1) does not depend 
on x.) The needed similarity is realized by the diagonal matrix 



A 



N 



= diag f(l - abq 2 ), 1(1 - abq A ), 1(1 - a6 9 6 ), . . . , -J-^l - abq 2N )) 
V 9 9 9 J 



so that 



AJ^SnMnS^An = -{Qn + A(s - 1)/jv) - A(x)7 



AT • 



(19) 



Verification of (|19|) is straightforward for off-diagonal entries. For diagonal entries, it reduces to verification 
of the identity 



cn+i - 1 - 



bq 



N 



1 

q + bq~ W 



9 = a n+ i 



1 — abq 



2n+4 



1 — a&q 



2n+2 



2 l-a6 9 2 



1 - a6 9 2 "+ 2 : 



where a n are given by and c„ by (|14|l . This last identity can be checked using M ATL AB Symbolic Math 
Toolbox. This finishes the proof of (JTSJl and |JTSJ . □ 
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The determinant for Racah polynomials is therefore 

RA N+1 (X(x); a, (3, 7, N) = (-x) N+1 (x + 7 + S + l) N +i 
and the determinant for Hahn polynomials is 

H N+1 (\(x);a,@, N) = (-x) N+1 . 
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